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Helium atom scattering (HAS) studies of the H-covered 
Mo (110) and W(110) surfaces reveal a twofold anomaly in 
the respective dispersion curves. In order to explain this un- 
usual behavior we performed density-functional theory calcu- 
lations of the atomic and electronic structure, the vibrational 
properties, and the electronic susceptibility of those surfaces. 
Our work provides evidence for hydrogen adsorption induced 
Fermi-surface nesting. The respective nesting vectors are in 
excellent agreement with the HAS data and recent angle re- 
solved photoemission experiments of the H-covered alloy sys- 
tem Moo.95Reo.05 (HO). Also, we investigated the electron- 
phonon coupling and discovered that the Rayleigh phonon 
frequency is lowered for those critical wave vectors compared 
to the clean surfaces. Moreover, the smaller indentation in 
the HAS spectra can be clearly identified as a Kohn anomaly. 
Based on our results for the susceptibility and the recently 
improved understanding of the He scattering mechanism we 
argue that the larger anomalous dip is due to electron-hole 
excitations by the He scattering. 



I. INTRODUCTION 

Much interest and many questions have been raised by 
the recent discovery of deep and extremely sharp inden- 
tations in the surface phonon spectra of H/Mo (110) and 
H/W(110) at full hydrogen coverage |||| (see Fig.0). 
Those anomalies are seen at an incommensurate wavevec- 
tor, Q ' , along the [001] direction (TH) and additionally 

at the commensurate wavevector Q c ' = S at the bound- 
ary of the surface Brillouin zone (SBZ). Out of the ordi- 
nary Rayleigh mode two simultaneous anomalies develop 
at those points. One, oj\, is extremely deep, and is only 
seen by helium atom scattering (HAS). The other, u>2, is 
instead soft, and is observed by both HAS and high reso- 
lution electron energy loss spectroscopy (HREELS) [§,@. 

A pronounced but less sharp softening of surface 
phonons is also seen on the clean (001) surfaces of W |j| 
and Mo ||. There, the effect can be explained in terms 
of a nested structure of the Fermi surface However, 
for the phonon anomalies of H/W(110) and H/Mo(110) 
such a connection between the electronic structure and 
the vibrational properties seemed to be non-existent since 
angular resolved photoemission (ARP) experiments of 
the two systems gave no evidence for nesting vectors 
comparable to the HAS determined critical wave vec- 
tors |]|l(| • An alternative model which links the phonon 
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FIG. 1. HREELS spectra of the clean W(110) surface (up- 
per panel) and its H-covered (1 x 1) phase (lower panel) plkj: 
Shown are the dispersion of the Rayleigh wave (triangles) and 
the longitudinally polarized surface phonons of the first (cir- 
cles) and second (squares) layer. The dots represent the re- 
sults of the HAS measurements Jll| 



anomalies to the motion of the hydrogen adatoms has 
to be ruled out because the HAS spectra remain practi- 
cally unchanged when deuterium is adsorbed instead of 
hydrogen |l|,p|. 

Another puzzle is added by the observation of a sym- 
metry loss in the low-energy electron diffraction (LEED) 
pattern of W (110) upon H-adsorption. The phenomena 
may be caused by a H-induced displacement of the to 
layer W atoms along the [110] direction jll) (see Figs 
and c). Similar studies for H/Mo(110) do not pro vide 
any evidence for a top-layer-shift reconstruction |12fl . 

Within a hrst-principles approach, the goal of our work 
is to give a comprehensive explanation for the observed 
anomalous behavior and clarify the confusing picture 
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drawn by the different experimental findings. 
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FIG. 2. (a) Adsorbate positions within the p(l x 1) unit 
cell of the W(110) and Mo (110) surfaces: Shown are the 
long-bridge (L), short-bridge (S), hollow (H), and on-top (T) 
site. (b-c) Structure of the H-covered W(110) surface 
as suggested by Chung et al. juj: (b) low H-coverage 
O < 0.5 ML; (c) top-layer-shift reconstruction for higher 
H-coverages Q > 0.5 ML. The white (shaded) circles repre- 
sent the W/Mo atoms of the surface (subsurface) layer. The 
H positions are indicated by full dots, (d) Structure parame- 
ters used in Table |. 



II. METHOD 

We perform density-functional theory (DFT) calcula- 
tions using the local-density approximation (LDA) for 
the exchange-correlation energy functional For 
the self-consistent solution of the Kohn-Sham equations 
a full-potential linearized augmented plane-wave (FP- 
LAPW) code jl4j is employed which we enhanced by 
the direct calculation of atomic forces |l5|,|l6|]. Within 
a damped molecular dynamics approach this enables an 
efficient determination of fully relaxed atomic geometries 
and frozen-phonon energies. For methodical details we 
refer to Ref |l7| ). In the case of the W surfaces the va- 
lence and semi-core (core) electrons are treated scalar 
(fully) relativistically while for Mo all electrons are han- 
dled non-relativistically. The in-plane lattice constants 
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= 3.14 A and 
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'Mo 



= 3.13 A are calculated without 



including zero-point vibrations. They are in good agree- 
ment with the respective measured bulk lattice parame- 
ters at room temperature which are 3.163A and 3.148A 
for W Ul and Mo @, respectively. 

We performed systematic tests comparing the 
LDA and the generalized gradient approximation |2C[] 
exchange-correlation functionals. For the quantities re- 
ported in this paper (total energy differences of different 
adsorption sites, bond lengths, etc.) both treatments give 
practically the same results. Also, we checked the kn- 
point convergence: In the case of atomic and electronic 
structure calculations a two-dimensional linear mesh of 
64 k|| -points within the (lxl) SBZ gives stable results. 
For the evaluation of frozen-phonon energies a set of 
56k||-points is employed within the SBZ of the enlarged 
(1 x 2) and (2 x 1) surface cells. 



III. RESULTS 
A. Atomic Structure 

The first step in our study is to determine the atomic 
structures of the clean and hydrogen-covered (110) sur- 
faces. The relaxation parameters presented in Table || 
are calculated for a nine-layer slab (the parameters are 
defined in Fig. ||d). The results, which are remarkably 
similar for Mo and W, are in excellent quantitative agree- 
ment with a recent LEED analysis [pl|p2[ . 

As the energetically most favorable hydrogen adsorp- 
tion site on both W(110) and Mo (110) we identify a 
quasi-threefold position (indicated as H in Fig.^Ja). Our 
investigations also throw light on the suggested model of 
a H-induced structural change: For both materials the 
calculated shift yi is only of the order of 0.01 A and thus 
there is no evidence for a pronounced top-layer-shift re- 
construction. Moreover, this subtle change in the surface 
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TABLE I. Relaxation parameters for the clean and 
H-covered (110) surfaces of Mo and W (see also Fig.^d). For 
each system the results from a nine-layer slab calculation (first 
line) and a LEED analysis (second line) |5l],^2| are presented. 
The height of the hydrogen above the surface and its [TlO] off- 
set from the [001] bridge position are denoted by du and j/h, 
respectively. The shift of the surface layer with respect to the 
substrate is j/i. The parameters Adij describe the percent- 
age change of the interlayer distance between the i-th and the 
j-th substrate layers with respect to the bulk interlayer spac- 
ing do- The numerical accuracy of the theoretically obtained 
parameters is ±0.015 A^ ±0.7%d - 
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2/i Adi2 

(A) (%d ) 
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Ad 34 
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Mo (110) 

-5.0 


±0.7 


-0.3 






-4.0 ± 0.6 ±0.2 ± 0.8 0.0 ± 1.1 


0.62 


1.09 


H/Mo(110) 
0.04 -3.3 


±0.3 


-0.2 


0.55 ±0.4 


1.3 ±0.3 


0.0 ±0.1 -2.0 ±0.4 


0.0 ±0.5 


0.0 ±0.8 






W(110) 

-3.6 


±0.2 


-0.3 






-3.1 ±0.6 


0.0 ±0.9 


0.0 ± 1.0 


0.67 


1.09 


H/W(110) 
0.02 -1.4 


-0.3 


-0.1 


0.56 ±0.4 


1.2 ±0.25 


0.0 ±0.1 -1.7 ±0.5 


0.0 ±0.6 


0.0±0.9 



geometry is unlikely to be resolved experimentally due 
to zero-point vibrations. This aspect becomes evident 
by evaluating the total energy with respect to a rigid 
top-layer-shift. In Fig. || we present the results of such 
a calculation for H/W(110) and depict the first three 
vibrational eigenstates obtained from a harmonic expan- 
sion of the total energy E tot (y 1 ); for H/Mo(110) the en- 
ergetics is similar. Since we have k^T ps 25meV at room 
temperature thermal fluctuations of y\ are of the order of 
0.1 A. This is considerably larger than the theoretically 
predicted absolute value of y\. 
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FIG. 3. Change of total energy E tot versus top-layer-shift 
i/i : For each data point the whole surface is re- 
laxed keeping only the substrate [110]-coordinates y\ and 
3/2 = J/3 = A fixed. Also shown are the first three oscil- 
lator eigenstates calculated from a harmonic expansion of the 
total energy E tot (yi). 



high density of states at the Fermi level. More impor- 
tant, the new Fermi contour shows pronounced nesting 
features (see Fig. ||c). The magnitudes of the calculated 
nesting vectors and HAS measured critical wavevectors 
along TS and TH are listed in Table ||. For both sys- 
tems the agreement between theory and experiment is 
excellent. 

Our results are at variance with the photoemission 
studies by Kevan and coworkers |9|,[l0| (see Fig. 2 in 
Rcf. p7[]). However, in view of more recent ARP mea- 
surements [24] which we compare to our results in Fig. ^ 
we dare to suggest that those differences may be due to 
problems in the experimental analysis. The later work 
deals with the (110) surface of the alloy Moo.95Reo.05 but 
the surface physics of both systems Moo.95Reo.05 (HO) 
and Mo (110) should be practically the same because in 



B. Electronic Structure 

The similarities between Mo (110) and W(110) found 
for the surface geometries continue when we compare the 
electronic structures. For both systems the H adsorption 
alters the surface potential and induces the shift of a 
band with (d^ z 2_ r i,d xz ) character to higher binding en- 
ergies |l^j2^]. Due to this effect which is illustrated in 
Fig. |J the Fermi line associated with this band is moved 
into the band gap of the surface projected band structure, 
and the respective states become true surface states. 

The shifted (d^ z 2_ r 2 1 d xz ) band is characterized by a 



TABLE II. Theoretical Fermi-surface nesting vectors com- 
pared to critical wavevectors obtained by HAS and HREELS 
experiments [Q^|3). 
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theory 


experiment 


TH 


H/Mo(110) 


0.86 


0.90 




H/W(110) 


0.96 


0.95 


TS 


H/Mo(110) 


1.23 


1.22 




H/W(110) 


1.22 


1.22 
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FIG. 4. Band structure (upper figure) and Fermi sur- 
face (lower figure) of clean Mo (110) (dotted curves) and 
H/Mo(110) (solid curves). The shaded areas represent the 
projection of the bulk Fermi surface onto the SBZ of the (110) 
surface. 
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FIG. 5. Theoretical Fermi surfaces of the (a) clean and (c) 
H-covered Mo (110) surface. The solid (dotted) lines denote 
surface resonances or surface states which are localized by 
more than 60 % (30 %) in the MTs of the two top Mo layers. 
Shaded areas represent the (110) projected theoretical Mo 

c i c 2 

bulk Fermi surface. The arrows Q ' and Q ' are the critical 
wave vectors predicted by the calculations. Also presented 
are data points (+) which stem from an ARP study of the 
(b) clean and (d) H-covered Moo.95Reo.05 (HO) surface |E4j. 



the HAS anomalies is now verified by experiment. 



both cases the top layer consists only of Mo atoms. This 
assumption is also backed by test calculations where we 
simulated a Moo.95Tco.05 alloy be the virtual crystal ap- 
proximation and found only minor quantitative changes. 

From Fig. || it becomes clear that in particular for 
the important (d3 Z 2_ r 2, d xz ) surface band, which was not 
seen by Kevan's group, experiment and DFT now agree 
very well. There are however differences: Theory pre- 
dicts bands centered at S which are not seen by ARP 
whereas the band circle at V in Fig. ||d is only observed 
experimentally. Also, in the calculations we find an ellip- 
tical band centered at the T point whereas in ARP the 
same band has the shape of a rectangle. Those discrep- 
ancies are probably due to the fact that it is rather dif- 
ficult in theory (and experiment!) to clearly distinguish 
between surface states and surface resonances. Addition- 
ally, due to interactions between the two surfaces within a 
slab system the k|| -space location of a calculated surface 
resonance which is not strongly localized at the surface 
cannot be determined accurately. Furthermore, we note 
that DFT is not expected to describe the Fermi surface 
accurately. Nevertheless, the important theoretical find- 
ing of a H-induced surface state which can be related to 



C. Vibrational Properties 

At this point, the following question arises: How does 
the surface react to the apparent electronic instability 
due to the strong nesting features? In the case of nesting 
the surface phonons take the burden to furnish energy 
and momentum to the scattering process of electrons at 
the Fermi level. This leads to a softening of the related 
phonons. If the electron-phonon coupling is strong and 
the energetic cost of a surface distortion is small this 
softening could even trigger a reconstruction combined 
with the build-up of a charge-density wave as in the case 
of the (001) surfaces of W || and Mo fl. It is clear that 
one needs to perform frozen-phonon calculations in order 
to determine the actual str eng th of the coupling and the 
induced phonon softening [j25|. 

Fortunately, at 0°' = S such a calculation is par- 
ticularly convenient, since at the zone boundary point S 
the Raylcigh-wavc polarization is purely vertical, and the 
second layer is immobile by symmetry. The calculated 
frequencies for the S'-point and a second zone boundary 
point, N, are presented in Table III. For the latter our 
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TABLE III. Comparison of calculated frozen-phonon en- 
ergies and experimental values obtained by HAS and 
HREELS |!|. The theoretical phonon energies are obtained 
using a five-layer slab. Their numerical accuracy is about 



±lmeV. 



phonon 


system 


E ph 


(meV) 






theory 


experiment 


N 


W(110) 


15.4 


14.5 




H/W(110) 


17.6 


17.0 


S 


Mo (110) 


22.7 


~21 




H/Mo(I10) 


17.2 


<16 




W(110) 


18.3 


16.1 




H/W(110) 


12.0 


11.0 



(a) original Brillouin zone 




(b) backfolding (TS' = |rS) 
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results reproduce the experimentally observed increase of 
the Rayleigh wave frequency as hydrogen is adsorbed. At 
the S point we find that those stiffening effects are over- 
compensated: The strong coupling to electronic states at 
the Fermi level leads to a lowering of the phonon energy 
in good agreement with the experimental results. 

In Fig. U we schematically illustrate the mechanism 
which is responsible for this effect. It is, in fact, a text- 
book example of a Kohn anomaly due to Fermi-surface 
nesting (2(],|27j] : (a) Within the unperturbated system the 
(d 3z 2_ r 2,dxz) band cuts the Fermi level exactly midway 
between T and S. (b) The S'-point phonon couples to 
the states at the Fermi level and causes the backfolding 
of the SBZ. (c) The nuclear distortion associated with the 
S'-point phonon modifies the surface potential and hence 
removes the degeneracy at the new zone boundary point 
S'. The occupied states are shifted to lower energies. 
This amounts to a negative contribution of the electronic 
band structure energy to the total energy and thus a 
lowering of the phonon energy. 

A frozen-phonon study for the second nesting vector 

along TH is not performed because Q ' is highly non- 
commensurate; a calculation would be very expensive. 
However, since the character of the (d^, z 2_ r 2 1 d xz ) band 
does not change when shifting from the TS to the TH 
nesting we expect similar results. 

The calculation of the electron-phonon interaction at 
the Mo (110) and W (110) surfaces and its change due to 
hydrogen adsorption pinpoints the phonon character of 
the small anomaly ll>2 and identifies the interplay between 
the electronic structure and the vibrational spectra of the 
transition metal surfaces. Thus, our results which were 
recently confirmed by Bungaro |2q| within the framework 
of perturbation DFT clearly support the interpretation 
that the small dip observed by both HAS and HREELS 
is due to a Kohn anomaly. 
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FIG. 6. Schematic representation of the mechanism which 
causes the Kohn anomaly on H/Mo(110) and H/W(110). 
Shown is the bandstructure along SYS parallel to the nesting 
vector Q c ' in Figs.[B|. The form of the (d 3z 2_ r 2,d xz ) band is 
indicated by dashed lines. 

D. Electron- hole excitations 

For the deep and sharp anomaly lo\ the interpretation 
is less straightforward, particularly since it is only seen 
in the HAS spectra. One needs to understand the nature 
of rare-gas scattering in order to explain this unusual 
behavior. 

In a recent study we found that those scattering pro- 
cesses are significantly more complicated (and more in- 
teresting) than hitherto assumed p9]: The reflection of 
the He atom happens in front of the surface at a dis- 
tance of 2-3 A. More important, it is not the total elec- 
tron density of the substrate surface which determines 
the interaction but the electronic wave functions close 
to the Fermi level. In the case of the H/W(110) and 
H/Mo(110) systems it is thus plausible to assume that 
the He atom couples directly to the {d iz 2_ r 2, d xz ) surface 
states mentioned above and excites electron-hole pairs. 

In order to analyze the spectrum of those excitations 
in some more detail we evaluate the imaginary part of the 
electronic susceptibility x(qy , hco) within random phase 
approximation (RPA): 

im x(qn,M = 
A f 

= Im JZ2 / dk ll W k ||+q|| W k|| (/ k||+q|t - / k|| ) X 
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£k||+q N - £k|| — fib) — 17] 

/k„)x (1) 



A f 

= JZ2 / dk ll ^n+qii^ii (/k ll+q 



with -A? 



■^2 being the ky -space density. In our approach 
we use all eigenvalues ek„ and occupation numbers /ky 
obtained via a nine-layer slab calculations. Also, we take 
into account the localization Wk., of the respective state 
at a distance of 2.5 A in front of the substrate surface. 

The results obtained for Mo (110) and H/Mo (110) are 
presented in Fig. [?]: For the clean surface we find that 
the intensity of the electron-hole excitations decreases 
continuously as we move away from the T point towards 
the zone boundaries. This relatively smooth behavior of 
the susceptibility is modified considerably when H is ad- 
sorbed on the clean surface: Pronounced peaks appear, in 
particular one which stretches from the TH line through 
the SBZ to the symmetry point S. This is in excellent 
agreement with the HAS results. Also, there is a peak 
located at the second half of the TN symmetry line. It 
cannot be associated with an anomaly in the HAS spec- 
trum. This result may be due to the fact that selection 
rules for the interaction between the scattering He atoms 
and the electron-hole excitations are not considered. 

In view of the simplicity of the approach, the calculated 
spectrum of electron-hole excitations however seems to 
be well-connected to the HAS measurements. Moreover, 
it is surprising how detailed the Fermi surface nesting 
is reflected in the susceptibility at such a large distance 
from the surface where the HAS scattering process takes 
place. This finding supports the earlier suggestion fHf] 
that the giant indentations, seen only in HAS, are due 
to the excitation of electron- hole pairs. By contrast, 
in HREELS experiments the electrons are directly scat- 
tered at the atoms and thus only measure the phonon 
spectrum. Therefore, the strong anomaly is invisible for 
HREELS. 



IV. SUMMARY 

In conclusion, the major results of our work are as fol- 
lows: We demonstrate that both the W(110) and the 
Mo (110) surface do not show a pronounced top-layer- 
shift reconstruction upon adsorption of hydrogen, a re- 
sult confirmed by a recent LEED analysis M]. Also, we 
give a consistent explanation for the H-induced anomalies 
in the HAS spectra of W (110) and Mo (110). The H ad- 
sorption induces surface states of (d 3z 2_ r 2 , d xz ) character 
that show pronounced Fermi-surface nesting. The soft 
anomalous dip is clearly identified as a Kohn anomaly 
due to those nesting features | p7| , p3|j2^ ] . For an under- 
standing of the huge dip one needs to take into account 
that the scattering of rare-gas atoms is crucially influ- 
enced by interactions with substrate surface wave func- 
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FIG. 7. Imaginary part of the electronic RPA suscepti- 
bility x(q||,M of Mo (110) and H/Mo(110) calculated for 
fojj = 0.5 mRy « 7meV. The dashed line represents a value 
of 3 a.u. while each full line denotes an additional increase by 
1 a.u. (with increasing line width). The positions of the HAS 
measured anomalies within the SBZ are indicated by dots. 
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tions at the Fermi level p9| . The He atom couples ef- 
ficiently to the H-induced surface states on H/W(110) 
and H/Mo(110). We therefore conclude that the deep 
branch of the anomaly is predominantly caused by a di- 
rect excitation of electron-hole pairs during the scattering 
process. 

We hope that our first-principles study stimulates some 
more experimental work: For instance, the liquid-like be- 
havior of the H adatoms observed in the HREELS mea- 
surements || is still not understood. One would also like 
to know whether the HAS and HREELS spectra of the 
Moo.95Reo.05 (HO) alloy system reveal H-induced anoma- 
lies as expected. Finally, we call for a new ARP study 
of the clean and H-covered Mo (110) and W(110) sur- 
face and think that scattering experiments with atoms 
or molecules like Ne or H2 might provide additional in- 
sight into the interesting behavior of these surfaces. 

We thank Erio Tosatti for stimulating discussions. 
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